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Improved Systematic of pp Elastic Scattering Data 

V. Uzhinsk}!!! and A. GaloyaiH 
CERN, Geneva, Switzerland 

Unified systematic of elastic scattering data (USESD) proposed by the au- 
thors ( ar Xiv:1111. 4984 [hep-ph]) is based on symmetrized 2-dimensional Fermi 
distribution for pp elastic scattering amplitude in the impact parameter repre- 
sentation. It allows to describe differential cross sections of the reactions up 
to \t\ - 1.75 (GeV/c)2. To extend it to higher \t\ values we consider a two 
coherent exponential parametrization of the cross sections and show that it 
cannot describe the cross sections at small \t\ at Piab > 10 GeV/c. We extract 
a description of high \t\ region from the parameterization and couple it with 
USESD. As a result, we obtained a good description of all pp elastic scattering 
data at Piab > 10 GeV/c. It can be easily used in Glauber Monte Carlo codes 
for calculations of nucleus-nucleus interaction properties. 

Introduction 

Recently, authors of the paper [1] remembered the old parameterization of pp elastic scattering am- 
plitude: 



and applied it for a fitting of the Totem Collaboration data |3| on elastic pp scattering at 



(1) 



7 TeV. 



The amplitude was proposed in 1973 by R.J.N. Phillips and V.D. Barger. Nearly at the same time, it 
was independently proposed and tested in papers [5J O [B] where it was applied for a fitting of pp elastic 
scattering. 

The authors of the paper |2] analyzed only pp experimental data at Piab =12, 14.2, 19.2, 24, 29.7 
GeV/c and at y/Spp = 53 GeV. No /NoF and the parameter's errors were given by them. So, a quality 
of the parameterization is unknown. It's predictive power is unclear. 

The authors of the papers [HOIB] fitted antiproton-proton scattering experimental data at Piab =1.11 
- 16 [B] GeV/c (twenty-eight sets of data). A fitting of eight selected data sets gave smallest errors. Three 
coherent exponentials was proposed and used in the paper [7] . There is not any reclamation to technical 
details of the fittings. We repeated the work and published the results in the paper where we proposed 
energy dependencies of the parameters. 

An advantage of Eq. [1] is that it can be easily applied in Glauber model calculations of hadron- 
nucleus and nucleus-nucleus interaction properties at high and super high energies, and it can improve 
an exactness of the calculations. Usually, only one exponent is used in the calculations. The aim of our 
present paper is a checking of Eq. [T]in a wide range of energies. 

There is another parameterization of high energy pp elastic scattering data proposed by us in the 
paper [9]: 



ndq Ji{Rq) 
■ih{Trdq) Rq 



1 ndq / T:dq 



2g2 sinh{'Kdq) \tanh{'Kdq) 



-1 MRq) 



q = V~t. (2) 



It is a Fourier-Bessel transform of a symmetrized 2-dimensional Fermi- function |10j . 



7(6) = A 



1 



1 



l_^g(6-fl)/d l^p-{b+R)/d 



- 1 



(3) 



where 7(6) is the elastic scattering amplitude in the impact parameter representation. 

A possibility to describe a high Pt elastic pp scattering was considered in Ref. [3] following papers 
HmH]; but we were not satisfied by results. Thus, below we undertake an effort to combine our approach 
with two coherent exponential one. As a result, we have obtained an improved parameterization of pp 
elastic scattering data. 
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1 Validation of two coherent exponential expression 

A direct application of Eq.[T]is complicated by strong parameter correlation. Thus, we write a differential 
elastic scattering cross section using Eq. [T] as: 



da 



(4) 



and fit experimental data at Piab > 10 (GeV/c). We took only 45 from 64 sets of experimental data 
which gave meaningful results. Fit results of the selected data sets are presented in Fig. 1, 2. 
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Figure 1: Energy dependencies of the parameters Ai and Bi. 
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Figure 2: Energy dependencies of the parameters A2, B2 and A3. 

As seen, the results are rather unstable. Though, x^/NoF =2572/1803~1.43. Of course, selecting 
only some experimental data one can obtain stable values of the parameters, as it was done in Ref. 

More serious drawback of the parameterization[T]is that it cannot describe small angle scattering data, 
especially, total cross sections, atot, and ratio of real and imaginary parts of elastic scattering amplitude 
at t = 0, p{0) = RefiO, s)//m/(0, s). 
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Figure 3: atot and Ref /Imf as functions of energy. Black points are experimental data from PDG data 
base [T3j. Squares are calculation results without error bars. 



Thus, we cannot recommend to use the parameterization for practical calculations. 
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2 Combination of the approaches 



Another situation takes place with the parameterization [51 It describe small angle scattering data rather 
well, but falls down at large \t\. At the same time, the parameterization [T] describes large \t\ data quite 
well. Thus, taking into account that (/> ^ tt in Eq. [T]we can combine the parameterizations as: 



Imf{s,t) = Ai 



^2 T^dq JijRq) ^ 1 T^dq f -ndq v Jo(i?g) 

sinhindq) Rq sinh{TTdq) \tanh{'!Tdq) 



(5) 



Ref{s,t)=A,-p-{Ry2 + n^d'/6)-^^^^JoiRq) + e'^^^'/^^ (6) 
p = 0.135 -3/Vs + 4/s + 80/s^ (7) 

^ = 10 • 25.68185 • TT • (/to/2 + i?e/2) \mb/(GeV/cy], (8) 
dt 

where q = V—t 25.68185 [fm^-^], t is 4-momentum transfer [{GeV/c)'^], R and d are in [/m], A2 is 
in [/to^]. The parameter Ai was introduced in order to take into account uncertainties in absolute 
normalization of experimental data. Experimental data on p [13] were approximated by Eq. [71 

There are a lot of experimental data, but most of them are small angle ones. They do not allow an 
unambiguous determination of the parameters. Thus, we select data at piab =14.2, 19.2, 24, 200, 293, 
501 GeV/c [m, and =44.7, 52.9, 62.5, and 7000 GeV [HI [^ in which regions of minimum at \t\ ~ 
1.5 (GeV/c)2 and high \t\ tails are presented. 

Results of the fit {x^ /NoF = 1.51) of the parameters A2 and B2 are presented in Fig. 4. 
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Figure 4: Points are fit results for A2 and i?2. Lines are approximations of the energy dependencies. 
The energy dependencies of the parameters can be approximated as: 

= 1.77 10"-^ [log(s/225)]2 + 0.05/s, B2 = 0.283 * log (s) + 30/s - 0.75. (9) 

A change of the sign of A2 in Eq. [Blfrom "+" to "-" makes the fit worse. 

Having the result for high \t\ values we can now determine "soft" parameters {Ai, R and d) more 
exactly than it was done before. In our paper ^ we found the parameters at an artificial restriction on 
"soft" amplitude apphcation region (Eq. [5]), \t\ < 1.75 (GeV/c)^. Now we fix A2 and B2 by Eq. [HI and 
fit other parameters using 64 sets of experimental data (see references in [9]). Results of the fitting are 
presented in Fig. 5 (x^/NoF =1.96). 

As seen, energy dependencies of R and d are determined rather well. They can be parameterized as: 

= 0.9/s°-2^ + 0.053 log(s), d = 0.379 - 0.26/s°-2^ (10) 

Fitted values of Ai scatter rather strong especially at low energies. We propose the following param- 
eterization of its energy dependence: 

Ai = 1.077 - 0.175 e-^ o^i"" ' + 0A5/s° '^^, (11) 
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According to Eq. [2] total cross section is given as: 

atot^2^ Ai {R" + Ti^(f/i). (12) 

If we use only fitting results for Ai, R, and d, we obtain green points in Fig. 6 shown without error bars. 
As seen, they are scattered rather strong. 




Figure 5: Points are fit results. Lines are approximations of the energy dependencies. 
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Figure 6: Open and closed points are experimental data from PDG [T? data base. Lines are calculation 
results using approximations of Ai(see text). 

Using the approximations given by Eqs. [71 [SI and [TUj {Ai was not fixed) we have /NoF — 
8179/2256 ~ 4.79. A fixing of Ai by Eq. [TT] leads to an increasing of ^ jNoF in 2 times. A qual- 
ity of the experimental data descriptions is presented in Appendix A. 



3 Application of USESD in Glauber Monte Carlo codes 

Glauber Monte Carlo codes [111 [T7] (see review in [TH]) calculate various properties of inelastic nucleus- 
nucleus interaction such as: number of participating nucleons, multiplicity of intra-nuclear collision, 
impact parameter distributions and so on. Mainly they use so-called "nucleon inelastic overlap function", 
(?(&), in the simplest form: 

g(h) = 0{Rint - b), Rint = ^/^i^- (13) 
In the Glauber theory the function is given as: 

g(6)=7(6)+7*W-|7(&)r (14) 

According to Eqs. [SI [H 

r ^-{h-R)/d I ^ 



4 



-i pAi- 



Rd 



1 _^g-(&-fl)/d 
A. 



2 



-e 



1 + g~{b+B^)/d 
b-'/{2 B2 25.64) 



2 



2ttB2 25.64 

= + (0.07 + ff + 0.2^2) e"^-^^/''. (16) 

A reason of the introduction of R is that Eq. [5] was obtained from Eq. [3] assuming d/R << 1. As 
seen from the fitting resuhs it is not true for pp interactions. Thus, the correction of Eq. was found 
at numerical investigation. 

Conclusion 

USESD is enlarged by the simple description of high |i| elastic scattering. The new parameters have been 
determined. Exact formulae are presented. 
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Appendix A: Comparison of experimental 
ization 




10"^ I , 1 , 1 , 1 , 1 , L 

1 2 3 4 5 

|t| (GeV^) 



Figure 7: The points are the experimental data 
by R.M. Edelstein et al., Phys. Rev. D5 (1972) 
1073. 
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Figure 9: The points are the experimental data 
by G.W. Brandenburg et al., Phys. Lett. 58B 
(1975) 367. 



a on pp-interactions with USESD parameter- 
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Figure 8: The points are the experimental data 
by J.V. Allaby et al., Nucl. Phys. B52 (1973) 
316. 
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Figure 10: The points are the experimental data 
by K.J. Foley et al., Phys. Rev. Lett. 11 (1963) 
425. 
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Figure 11: The points are the experimental data 
by K.J. Foley et al., Phys. Rev. Lett. 15 (1965) 
45. 
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Figure 13: The points are the experimental data 
by D. Harting, Nuovo Cimento 38 (1965) 60. 
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Figure 12: The points are the experimental data 
by J.V. Allaby et al., Nucl. Phys. B52 (1973) 
316. 
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Figure 14: The points are the experimental data 
by K.J. Foley et al., Phys. Rev. Lett. 11 (1963) 
425. 
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Figure 15: The points are the experimental data Figure 16: The points are the experimental data 
by J.V. Allaby et al., Nucl. Phys. B52 (1973) by K.J. Foley et al, Phys. Rev. Lett. 11 (1963) 
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Figure 17: The points are the experimental data 
by R.M. Edelstein et al, Phys. Rev. D5 (1972) 
1073. 
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Figure 18: The points are the experimental data 
by .K.J. Foley et al, Phys. Rev. Lett. 11 (1963) 
425 



P, =19.2 GeV/c 

lab 




|t| (GeV^) 

Figure 19: The points are the experimental data 
by D. Harting, Nuovo Cimento 38 (1965) 60. 
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Figure 20: The points are the experimental data 
by J.V. Allaby et al., Phys. Lett. B28 (1968) 67. 
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Figure 21: The points are the experimental data 
by .K.J. Foley et al., Phys. Rev. Lett. 11 (1963) 
425 
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Figure 22: The points are the experimental data 
by R.M. Edelstein et al., Phys. Rev. D5 (1972) 
1073. 
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Figure 23: The points are the experimental data 
by J.V. Allaby et al., Phys. Lett. B28 (1968) 67. 
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Figure 25: The points are the experimental data 
by R.M. Edelstein et al, Phys. Rev. D5 (1972) 
1073. 
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Figure 27: The points are the experimental data 
by C. Bruneton et al., Nucl. Phys. B124 (1977) 
391. 
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Figure 24: The points are the experimental data 
by J.V. Allaby et al., Nucl. Phys. B52 (1973) 
316. 
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Figure 26: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 28: The points are the experimental data 
by Z. Asad et al., Nucl. Phys. B255 (1984) 273; 

C. W. Akerlof et al., Phys. Rev. D14 (1976) 2864; 

D. S. Ayres et al., Phys. Rev. D15 (1977) 3105. 
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Figure 29: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 30: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 31: The points are the experimental data 
by D.S. Ayres et al., Phys. Rev. D15 (1977) 
3105. 



100 GeV/c 




|t| (GeV^) 



Figure 32: The points are the experimental data 
by C.W. Akerlof et al., Phys. Rev. D14 (1976) 
2864; R. Rubinstein et al., Phys. Rev.DSO (1984) 
1413. 
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Figure 33: The points are the experimental data 
by D.S. Ayres et al., Phys. Rev. D15 (1977) 
3105; R. Rusack et al., Phys. Rev. Lett. 41 
(1978) 1632. 
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Figure 34: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
101632. 
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Figure 35: The points are the experimental data 
by D.S. Ayres et al., Phys. Rev. D15 (1977) 
3105. 
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Figure 36: The points are the experimental data 
by D. Brick et al., Phys. Rev. D25 (1982) 2794. 
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Figure 37: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 38: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 39: The points are the experimental data 
by D.S. Ayres et al., Phys. Rev. D15 (1977) 
3105. 
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Figure 40: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 41: The points are the experimental data 
by C.W. Akerlof et al., Phys. Rev. D14 (1976) 
2864; R. Rubinstein et al., Phys. Rev.DSO (1984) 
1413. 
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Figure 42: The points are the experimental data 
by A. Schiz et al., Phys. Rev. D24 (1981) 26; G. 
Fidecaro et al., Nucl. Phys. B173 (1980) 513. 
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Figure 43: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 45: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 



Figure 44: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 46: The points are the experimental data 
by M.G. Albrow et al., Nucl. Phys. B108 (1976) 
1. 
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Figure 47: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 
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Figure 49: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 
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Figure 51: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 
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Figure 48: The points are the experimental data 
by M.G. Albrow et al., Nucl. Phys. B108 (1976) 
1. 
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Figure 50: The points are the experimental data 
by A. Breakstone et al. Nucl. Phys. B248 (1984) 
253. 
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Figure 52: The points are the experimental data 
by E. Nagy et al., Nucl. Phys. B150 (1979) 221. 



13 



Vs= 53 GeV 




|t| (GeV^) 



|t| (GeV^) 



Figure 53: The points are the experimental data 
by A. Breakstone et al. Nucl. Phys. B248 (1984) 
253; Phys. Rev. Lett. 54 (1985) 2180. 
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Figure 54: The points are the experimental data 
by S. Erhan et al., Phys. Lett. B152 (1985) 131; 
J.C.M. Armitageetal.,Nucl. Phys. B132 (1978) 
365. 
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Figure 57: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 



Figure 58: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 
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Figure 59: The points are the experimental data 
by the Totem Collaboration (G. Antchev et al., 
Europhys. Lett., 96, 21002 (2011); Europhys. 
Lett., 95, 41001 (2011)) digitized by us. 
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Figure 60: The points are the experimental data 
by the Totem Collaboration (G. Antchev et al., 
Europhys. Lett., 96, 21002 (2011); Europhys. 
Lett., 95, 41001 (2011)) digitized by us. 
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Figure 61: The points are the experimental data by the Totem Collaboration (G. Antchev et al., Europhys. 
Lett., 96, 21002 (2011); Europhys. Lett., 95, 41001 (2011)) digitized by us. 
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